Measurement of the Light Antiquark Flavor Asymmetry in the Nucleon Sea 
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A precise measurement of the ratio of Drell-Yan yields from an 800 GeV/c proton beam incident 
on hydrogen and deuterium targets is reported. Over 140,000 Drell-Yan muon pairs with dimuon 
mass M M + M - > 4.5 GeV/c 2 were recorded. From these data, the ratio of anti-down (d) to anti-up 
(m) quark distributions in the proton sea is determined over a wide range in Bjorken-a;. A strong x 
dependence is observed in the ratio d/u, showing substantial enhancement of d with respect to u for 
x < 0.2. This result is in fair agreement with recent parton distribution parameterizations of the 
sea. For x > 0.2, the observed d/u ratio is much nearer unity than given by the parameterizations. 
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No known symmetry requires equality of the d and u 
distributions in the proton. Until recently it had been 
generally assumed that d(x) = u(x) for lack of experi- 
mental evidence to the contrary. This assumption may 
be evaluated by use of the expression 



[dp fx) - Up(x)] dx = 0.147 ± 0.039, 



(2) 



a considerable excess of d p relative to u p 



mx)-F-{ X )]^ = \- 2 - 



[d p (x) — u p (x)] dx. 

(1) 



Here F%(x) and F£{x) are the proton and neutron in- 
clastic structure functions, and d p (x) and u p {x) are the 
anti-down and anti-up quark distributions in the proton 
sea as a function of Bjorken-a;. Equation [l] requires the 
assumption of charge symmetry between the proton and 
neutron (i.e. u p = d n , u p = d n , etc.). If the nucleon sea is 
flavor symmetric in the light quarks, the value of the inte- 
gral on the left is 1/3, a result referred to as the Gottfried 
Sum Rule (GSR) |J. In 1991 the New Muon Collabora- 
tion (NMC) at CERN presented evidence that the GSR 
is violated, based on deep inelastic muon scattering data 
from hydrogen (p) and deuterium (d). They reported a 
final value of ^ [if (x) - F£(x)] = 0.235 ± 0.026 [§, 
which implies that 



This result has 
been adopted in the most current parameterizations of 
the parton distributions in the nucleon (^^]. 

Following publication of the NMC result, the use of 
the Drell-Yan process || was suggested || as a means by 
which the light antiquark content of the proton could be 
more directly probed. This was first done by the Fermilab 
E772 collaboration. They compared the production of 
Drell-Yan muon pairs from isoscalar targets to that from 
a neutron rich target. This measurement sets constraints 
on the non-equality of u and d in the range 0.04 < x < 
0.27 0. Later, the CERN experiment NA51 || carried 
out a comparison of the Drell-Yan muon pair yield from 
hydrogen and deuterium at a single value of x using a 
450 GeV/c proton beam and found 



u„ 



= 0.51 ±0.04 ±0.05. 



(3) 



(i;)=0.18 



A recent review by Kumano || presents an extensive dis- 
cussion of the existing literature on the flavor asymmetry 
of the antiquark sea. 
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Fermilab Experiment 866 (E866) measured the Drell- 
Yan muon pair yield from 800 GeV/c proton bombard- 
ment of liquid deuterium and hydrogen targets. From 
these data, d/u and d — u in the proton over the range 
0.020 < x < 0.345 are extracted. A significant difference 
between the d and u distributions is found. 

E866 used a 3-dipole magnet spectrometer em- 
ployed in previous experiments (E605, E772, and E789), 
modified by the addition of new drift chambers and ho- 
doscopes with larger acceptance at the first tracking sta- 
tion. Other improvements to the spectrometer included 
a programmable trigger system [jflj and a VME-based 
data acquisition system. An 800 GeV/c extracted pro- 
ton beam with up to 2 x 10 12 protons per 20 s spill 
bombarded one of three identical 50.8 cm long cylindri- 
cal target flasks containing cither liquid hydrogen, liq- 
uid deuterium or vacuum. After passing through the 
target, the remaining beam was intercepted by a cop- 
per beam dump located inside the second dipole magnet. 
The beam dump was followed by a 13.4 interaction length 
absorber wall of copper, carbon and polyethylene which 
blocked the entire aperture of the magnet. This absorber 
wall removed hadrons produced in the target and the 
dump, ensuring that only muons traversed the spectrom- 
eter's detectors. The detection system consisted of four 
tracking stations and a momentum analyzing magnet. 
The spectrometer's acceptance as a function of pr, the 
transverse momentum of the dimuon pair, was reasonable 
to 3.0 GeV/c, with some acceptance to 5.0 GeV/c. 

The targets alternated between hydrogen and deu- 
terium every five beam spills with a single spill collected 
on the empty flask at each target change. The tar- 
gets were 3.54 and 8.14 g/cm 2 thick, corresponding to 
7 and 15% of an interaction length for hydrogen and 
deuterium, respectively. Beam intensity was monitored 
by secondary-emission detectors, an ion chamber and 
quarter-wave RF cavities. Two four-element scintillator 
telescopes viewing the targets at nearly 90° monitored 
the luminosity, beam duty factor and data acquisition 
live time. The trigger required a pair of triple hodoscope 
coincidences having the topology of a muon pair from the 
target. Typically 70 triggers per second were recorded 
with an electronic live time above 98%. An integrated 
flux of 1.3 x 10 17 protons was delivered. 

Over 330,000 Drell-Yan events were recorded, using 
three different spectrometer settings which focused low, 
intermediate and high mass muon pairs. The data col- 
lected with the low and intermediate mass settings have 
systematic effects of a few percent which require addi- 
tional study. The data from the high mass setting are 
relatively free from these effects due to the greatly re- 
duced rates in the tracking chambers. Therefore, this 
Letter presents only the results from the high mass set- 
ting, with over 140,000 Drell-Yan events. 

In calculating the Drell-Yan yields, a small correction 
(averaging 0.2%) for random coincidences between two 
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FIG. 1. The ratio a pd /2a pp of Drell-Yan cross sections vs. 
X2- Shown are next-to- leading order calculations, weighted 
by acceptance, of the Drell-Yan cross section ratio using the 
CTEQ4M (solid) and MRS(R2) (dashed) parton distribu- 
tions. Also shown is a leading order calculation of the cross 
section ratio using CTEQ4M (dotted). In the lower CTEQ4M 
curve d — u has been set to as described in the text. The 
errors are statistical only. An additional 1% systematic un- 
certainty is common to all points. 

unrelated, oppositely charged muons was made. This 
correction was evaluated by studying the observed rates 
of same charge muon pairs. A background rate of approx- 
imately 8% from the target flask and beam line windows 
was measured using the evacuated flask and was sub- 
tracted. Given the identical geometry of the deuterium 
and hydrogen targets and the slight difference in the av- 
erage interaction point due to beam attenuation in the 
target, acceptance differences between the two targets 
were shown by a Monte Carlo simulation of the spec- 
trometer to be very small. The measured cross section 
ratio was corrected for differences in beam attenuation 
in the targets, target density and a small hydrogen con- 
tamination in the deuterium target. Additionally, the 
high luminosity caused a small, rate dependent ineffi- 
ciency. This led to a correction of 1.2% in the cross sec- 
tion ratio. The systematic error in the ratio of yields 
from the two targets is dominated by the uncertainties 
in the rate dependence (±0.6%), hydrogen contamination 
in the deuterium target (±0.2%) and beam attenuation 
(±0.2%). All other contributions, including the produc- 
tion of muon pairs from secondary hadron reinteraction, 
are small. The total systematic error in the cross section 
ratio is less than ±1%. 

The resulting ratio of the Drell-Yan cross section per 
nucleon for p + d to that for p + p is shown in Fig. |l| and 
Table | as a function of X2 [ p^JT^ I , the momentum frac- 
tion (Bjorken-x) of the target quark in the parton model. 
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TABLE I. Ratio of deuterium to hydrogen cross sections 
per nucleon vs. X2- The average of the kinematic variables 
for each bin is also tabulated. Note that (2:1) = (xf) + (xz). 
The errors are statistical only. An additional 1% systematic 
uncertainty is common to all points. 
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denoted by x\.) 



Muon pairs with mass, M^^- , below 4.5 GeV/c 2 or be- 
tween 9.0 and 10.7 GeV/c 2 were removed to eliminate 
contributions from the J/yb and T resonance families. 
The data clearly show that the Drell-Yan cross section 
per nucleon for p + d exceeds p + p over an apprecia- 
ble range in x 2 . Figure [l] also shows the predictions for 
next-to-leading order calculations [[uj of the cross section 
ratio, weighted by the E866 spectrometer's acceptance, 
using the CTEQ4M f| and MRS(R2) § parton distri- 
butions. The lower curve shows the predicted ratio for a 
modified CTEQ4M parton distribution which maintains 
the parameterization for d p +u p but sets d p — u p = 0. The 
data are in reasonable agreement with the unmodified 
CTEQ4M and the MRS(R2) predictions for x 2 < 0.15. 
It is clear that d p ^ u p in this range. Above x 2 = 0.15 
the data lie well below both parameterizations. 

The acceptance of the spectrometer was largest for 
x-p = x\ — x 2 > 0. In this kinematic regime the Drell-Yan 
cross section is dominated by the annihilation of a beam 
quark with a target antiquark. This fact, coupled with 
the assumption of charge symmetry between the neutron 
and proton and the assumption that the deuteron parton 
distributions can be expressed as the sum of the proton 
and neutron distributions, yields a simple approximate 
form of the Drell-Yan cross section ratio, 



pd 



2<jPP 



1 



I 1 -I- 1 dl d-2 \ 



1+*L 

u 2 



(4) 



The subscripts 1 and 2 denote that the parton distribu- 
tions in the proton as functions of x\ and x 2 , respectively. 
In the case that d = u, the ratio is 1. This equation il- 
lustrates the sensitivity of the Drell-Yan measurement to 
d/u and is valid only for x\ 3> x 2 . It does, however, 
imply an excess of d with respect to u for the data. Esti- 
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FIG. 2. The ratio of d/u in the proton as a function of x 
extracted from the Fermilab E866 cross section ratio. The 
curve is from the CTEQ4M parton distributions. The error 
bars indicate statistical errors only. An additional systematic 
uncertainty of ±0.032 is not shown. The result from NA51 is 
also plotted as an open box. 



mates of the nuclear effects in deuterium are significantly 
less than statistical errors shown in Fig. [l] Jl5[ . 

Some of the data, especially at higher x 2 , do not sat- 
isfy the x\ ^> x 2 criterion of Eq. |^. Consequently, d/u 
was extracted iteratively by calculating the leading order 
Drell-Yan cross section ratio using a set of parton dis- 
tribution functions as input and adjusting d/u until the 
calculated cross section ratio agreed with the measured 
value. In this procedure, the values for the d + u, va- 
lence and heavy quark distributions given by the global 
fits [e.g., CTEQ4M and MRS(R2)] were assumed to_be 
correct. In the beam proton, when x\ < 0.345, the d/u 
distribution was assumed to be the same as in the target 
proton. For x\ > 0.345, a constant value of 1 for d/u was 
used. Varying the high-xi value of d/u produced almost 
no change in the low x 2 bins, and less than a 3% change 
in the highest x 2 bin. This procedure was followed us- 
ing both the CTEQ4M and MRS(R2) parameterizations 
and negligible differences were seen. The extracted d/u 
ratio is shown in Fig. |^ along with the prediction made 
by CTEQ4M. A qualitative feature of the data, not seen 
in either parameterization, is the rapid decrease towards 
unity of the d/u ratio beyond x — 0.2. At x — 0.18, the 
extracted d/u ratio is somewhat smaller than the value 
obtained by NA51. Although the average value of Q 2 
(M 2 + ^_) is different for the two data sets, the change 
in d/u predicted by the parton distributions due to Q 2 
evolution is small. 

To address the GSR violation observed by NMC, the 
extracted d/u ratio is used together with the CTEQ4M 
value of d + u to obtain d — u. (Nearly identical results 
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FIG. 

(■0.345 



3. Fermilab E866 values for (a) d — u and (b) 
f U Ml> (d— u) dx' in the proton versus x. The curves rep- 
resent the corresponding values obtained by the CTEQ4M 
(solid) and MRS(R2) (dashed) parameterizations. The bar 
at 0.147 ± 0.039 on the left axis in (b) shows the result ob- 
tained by NMC for the integral from to 1. 



are obtained if MRS(R2) is used instead of CTEQ4M.) 
Since each bin in a; has a different average Q 2 , d — u was 
scaled to a common Q value of 7.35 GeV/c, the average 
for the entire data set. Based on this, the integral of 
d — u between a; mm and 0.345 is calculated. Both d — u 
and J min (d — u) dx are shown in Fig. ^. The integral 
reaches a value of 0.068 ± 0.007(stat) ± 0.008(syst) at 
x min = 0.02. This may be compared with the CTEQ4M 
and MRS(R2) parameterizations which have values of 
0.076 and 0.100, respectively, for the integral over the 
same region. Over the range 10~ 4 < x < 1, CTEQ4M 
gives a value of 0.108 for the integral, and MRS(R2) gives 
0.160. Above x — 0.345, it is unlikely there are signif- 
icant contributions to the d — u integral since the sea 
is relatively small in this region. Both CTEQ4M and 
MRS(R2) find this region contributes less than 0.002 to 
the total integral. It is clear, however, that significant 
contributions to the integral arise in the unmeasured re- 
gion below x = 0.02. 

Such a large dj u asymmetry cannot arise from pertur- 
bative effects [|16|. Most parameterizations of the parton 
distribution functions (e.g. CTEQ, MRS) simply assume 
a shape for d—u that accommodates the NMC and NA51 
results. It has been suggested 0,|l8| that including the 
effects of virtual mesons can account for the observed 



asymmetry and this appears |19-22| to be at least qual- 
itatively correct. 

In summary, this Letter reports a measurement of the 
Drell-Yan cross section ratio per nucleon of p + d to p+p. 
From this measurement the asymmetry of the light quark 
sea in the proton is extracted as a function of x. A feature 



of the present result is the reduction in d/u for x > 0.2. 
The current data are in qualitative agreement with NA51, 
but with a smaller d/u ratio at the single value of x which 
they measured. Over the range of x covered in this exper- 
iment (0.020 < x < 0.345), the integral of d — u reaches 
a value of 0.068 ± 0.007(stat) ± 0.008(syst). However, 
contributions to the integral arise from the unmeasured 
region below x = 0.020. 
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